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Abstract 

Two quantum quartic anharmonic many-body oscillators are introduced. One of them 
is the celebrated Calogero model (rational An model) modified by quartic anharmonic two- 
body interactions which support the same symmetry as the Calogero model. Another model 
is the three-body Wolfes model (rational G2 model) with quartic anharmonic interaction 
added which has the same symmetry as the Wolfes model. Both models are studied in the 
framework of algebraic perturbation theory and by the variational method. 
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This work is dedicated to the memory of Ian Kogan who died so young that is hard for 
me to imagine that he is not with us anymore. I knew him for about 30 years since the 
time when he appeared at ITEP Theory Division as a young, very briUiant student. Then 
for many years we were sitting in the next door offices at ITEP. Sometimes we talked on 
science being both intrigued by the transition from Quantum Mechanics to Quantum Field 
Theory. I am sure that Ian would be pleased to read the present article. 



I. INTRODUCTION 



Anharmonic oscillators play a crucially important role in contemporary physics since 
they model intrinsic anharmonic effects of the real world. The goal of the present work is to 
introduce and then to study a special type of quantum anharmonic oscillator - many-body 
anharmonic oscillators. One of them can be considered as an anharmonic perturbation of 



the celebrated many-body Calogero model (see or, in other words, the rational An model. 
Another is an anharmonic perturbation of the three body Wolfes model or, equivalently, the 
rational G2 model (see 1^). The first system describes n interacting particles on a line with 
fixed ordering with pairwise interaction, while the second one corresponds to three identical 
interacting particles on a line with fixed ordering with two- and three-body interactions. 
It is rather natural to impose a requirement that these anharmonic systems should possess 
the same symmetry properties as the original Calogero or Wolfes models: (i) translation 
invariance, (ii) permutation invariance, (iii) reflection symmetry with respect to a change of 
the sign of all coordinates. 

The one-dimensional anharmonic oscillator 

n = -— + x^ + Xx\ xeR, (1) 

is perhaps one of the most celebrated and the most studied problems in quantum mechanics. 
A systematic study was carried out by Bender- Wu in 1969-1973 in their seminal papers 0|. 
Even this simplest anharmonic oscillator possesses exceptionally rich properties: 

• divergent perturbation theory (PT) in A 0,01, i.e. 

E = OfcA'^ , Ofc oc ! , 
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• highly non-trivial but convergent PT in 1/A (strong coupling expansion) [3|, |5|, i.e. 

• analytic structure in A; all even (odd) eigenvalues are analytically related through 
square-root branch points, which accumulate to A = (see Fig.l) jsf. 

Hence by studying one eigenstate the whole family of eigenstates is explored! 



• 
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FIG. 1: Structure of singularities in the inverse coupling constant = X^^^^ on the first sheet 
of the Riemann surface of the ground state energy. Bullets denote the square-root branch points 
connected by cuts (vertical lines) 

So far very little (almost nothing) is known about eigenf unctions as functions of A. 
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II. ANHARMONIC CALOGERO MODEL 



The Hamiltonian of the A„_i anharmonic oscillator takes the form 

1=1 i>j •' i>j i>j 

where g > — |, u is the frequency, A > is the coupling constant with a factor (n + 6) which 
is introduced for a convenience and n = 2,3, . . .. This Hamiltonian describes a system of 
n identical particles situated on the straight line with pairwise interaction separated from 
each other by impenetrable barriers. The configuration space is 

— OO < Xi < X2 < ■ ■ ■ < Xn < OO . (3) 

If the coupling constant A = 0, the Hamiltonian Q corresponds to the celebrated Calogero 
model and the domain Q is nothing but the Weyl chamber of the root system. The 
ground state eigenfunction of the Calogero model is 

^f^(x) = A''{x)e-^^^ , g = iy{u-l) , (4) 

where A(x) = Y[i<i I 

the Vandermonde determinant and X2 = Yli>j{^i ~ ^jYi 
when the ground state energy is E^'' = uj{l + vn). 

In order to deal with translation invariance of many-body systems we replace the Carte- 
sian coordinates by the center-of-mass coordinate, Y = ,'Y^=i^j-> ^^e translation- 
invariant relative coordinates - the Perelomov coordinates |6| , 

?/i = Xi - - y , i = 1,2, . . . ,n , (5) 

n 

which obey the constraint X]j'=i Vj — 0) where Xi are the Cartesian coordinates. The coordi- 
nates (0) make sense as translation-invariant relative coordinates which measure a distance 
from the center of mass to a particle position. Since we consider a system of identical 
particles, permutation symmetry holds. In order to make manifest the permutation symme- 
try we introduce permutationally symmetric coordinates. The most convenient candidate 
is the invariants of the symmetric group. Eventually, we arrive at elementary symmetric 
polynomials of the arguments y (see Eq.dS}) as new coordinates 

{Xi,X2,...Xn) {Y,Tk{x) =Ok{y{x))\ fc=2,3,...n) . (6) 
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Here, 

il<«2<-<*fc 

are elementary symmetric polynomials. As an illustration let us present in explicit form the 
r coordinates for n = 2,3, putting — = Vi + 1/2 + ■ ■ ■ + Vn-i, 

• n = 2 T2 = -yl , 



• n = 3 T2 = -yl- ym -yl , = -ymiyi + 1/2) • 

It is easy to recognize that the r coordinates are nothing but a particular form of the 
Weyl invariant polynomials of the lowest degrees in the An-i root space. 
It can easily be shown that 

n 

^(xi - XjY = -2nT2 , (7) 

i>j 
n 

^{Xi-Xjf = 2{n + Q)rl-4:nri. (8) 

i>j 

These relations reveal a remarkable feature of the r coordinates - although the left- 
hand-side depends on all coordinates, the right-hand-side depends on a finite number 
of r's. Making a gauge rotation of the Hamiltonian ((21) with the Calogero ground state 
eigenfunction (j3)) as the gauge factor and re- writing the result in the r coordinates (jHl), we 
arrive at a strikingly simple expression after separating out the center-of-mass coordinate 
Y, 

hA = 2{^i'Y' (^A - 4'^) = /^Cal + XV, (9) 



1,3=2 ■' 1=2 



2 2n 

^2 ^4 



where 



(n-i + l)(l-j) . ., 

Z>max(lJ— j) 

Bi = ( — \- uj {n — i + 2){n — i + 1) rj_2 + 2uiTi , 



Explicit formulae for the first few coefficient functions are 



A22 = 2x2, B2 = 4:UT2 + {l + un){n-l) , 
A23 = 3r3, B3 = 6UJT2, 

A2i = 4r4, = 8ljt4^-\ — (1 + h'n)(n — 2)(n — 3) T2 , 

n 

A33 = 4r4-2(l-^)r2, ^34 = 5 Ts - 2 (^1 - ^) TaTg 
Am = Qtq + 2T2T4- 3(^1 - . 



In [3] it was demonstrated that at A = the gauge-rotated Calogero Hamiltonian hca\ 
(jni) has infinitely many finite-dimensional invariant subspaces 

Vk = {rrrr...T/"\0<Ep,<k), k = 0,l,2,.... (11) 

These spaces can be embedded one into another, 

P0CV1CP2C ...CPkC ... , 

thus forming an infinite flag (filtration) V . Hence one can say that the operator hcai preserves 
the flag V. Another property of /icai is the existence of a hidden gln-i algebra. The 
Hamiltonian /icai can be written as a second degree polynomial in generators of the gin 
algebra in the totally symmetric representation [k, 0, 0, . . . , 0), 

Ji~ = i = 2,3...n , 

OTi 
^ 

J^ij — '^i'o^^ i, J = 2, 3 . . . n , 



1=2 * 



Ji^ = r,:r° = x, (^r,^-A;), 2 = 2,3...n, (12) 

j=2 '^^ 

in such a way that the generators J7^^ do not appear. It is worth mentioning that for 
integer k, r? the generators (fT^ possess a common finite-dimensional invariant subspace 
Pfc. This is nothing but a finite-dimensional irreducible representation space of the algebra 
gl\^ in the realization (jl2j) . Therefore the flag V is made out of irreducible finite-dimensional 
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representation spaces of the algebra gin taken in the reahzation p2|) . It is worth emphasizing 
that the perturbation potential is itself an element of the representation spaces (|TT|). 

VpeVk, A; = 2, 3, ... . (13) 



A. Perturbation theory (generalities) 

Now let us consider the spectral problem for the operator h^, 

hx<p = 2e<j). (14) 
The spectral parameter e is related to the energies E' of Q by 

E = ) + e , 
and is related to the eigenfunction of (j21) through 

^(a;) = 0(x)^[,'^(x) , 

where '^^q\x) is the ground state eigenfunction of the Calogero Hamiltonian (j21) at A = 0. 
We develop perturbation theory for the equation (fT^ in powers of A, 

oo oo 

0=^0fcA\ e = ^e,A^ (15) 

fe=0 fc=0 PI 

which is in fact the Dalgarno- Lewis form of perturbation theory It is easy to derive an 
equation to find the fcth correction 

k 

(/iCai - 2eo)</)fc = 2 ^ ei </)fc_i - Vp (f)k-i . (16) 
1=1 

Following the theorem from ^9], as a consequence of the property (fT^ the perturbation 
theory ()15p is algebraic - all perturbation corrections (pk are polynomials in r's of finite 
degree. Hence the construction of perturbation theory is a linear algebraic procedure. 

B. Perturbation theory (concrete results) 

1. Ground state 

The ground state of the gauge-rotated Calogero Hamiltonian hcni is (see ©) 

00 = 1 , eo = . 
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(17) 



A simple analysis of equation ()16p shows that the first eigenfunction correction 0i G P2, 
hence it should be a second degree polynomial in the r's. After substitution of such an 
Ansatz into (fTH|) simple calculations give 

n{n — 1)(1 + z/n)[6 — z/(6 — 5n)] 



4(n + 6)cj2 



n 1 2 , ri[6 - z/(6 - Sri)] 

■ ^4 - 1— ^2 77 — , o — ^2 • (loj 



2{n + Q)uj Au ' 4(n + 6V2 
It is worth mentioning that the correction 01 depends on two r variables only, t^^/^. 

A similar analysis of equation (fTBj) shows that the second eigenfunction correction 02 ^ 
P4, hence it should be a fourth degree polynomial in the r's. After substitution of such an 
Anzatz into simple calculations give for the second correction 

^2 = \ai i\^2 5 [I50n + 36 + z/(5n - 6) (49n + 6) (19) 

.2„/i ni „2 



(n + 6)^ 



+ u^n{101n' - 245n+ 150 
2 '^('^ + 6) 2 



n(7n + 8 - Qun + bvn^) n[l^n + 6 + un{14n - 19)] 
''"4T2 TT-A ^4 



4cu3 8uj^ 
n(n-3) 2 (^^ + 6)2 ^ (n + 6)[4(5n + 3) + 3i/ri(5n - 6)] 3 

55n2 + 120n + 36 + vn{7An^ + 5ra - 114) + v'^n^{bn - Qf ^ 
+ 16Z7 

_ n[150n + 36 + z/(49n + 6) (5n - 6) + i/2n(101n2 - 245n + 150)] 

16^1^5 

It is worth mentioning that the correction 02 depends on four r variables only, r2,3.4,6- 

A similar analysis of equation (fT?)|l shows that the third eigenfunction correction 03 G Pe, 
hence it should be a sixth degree polynomial in the r's. After substitution of such an Anzatz 
into ()16|) after simple straightforward calculations we get for the third energy correction 

.3 [18 (36 + 144n + 215,,^) (20) 

32(n + 6)'^u;* 

-9z/ (72 + 504n + 772n=^ - 1033^^) 



+nz/2(2592 + 6948n - 16524/1^- 7529/1^) 
-nV^(3870 - 9297n + 752971^- 2052n^ 



It is worth mentioning that the correction 03 depends on six r variables only, r2,3,4,5,6,) 



It can be shown that the kth correction to the eigenfunction 0^ G V2k for 2k + 2 < n; 
hence it should be a 2A;th degree polynomial in the r's. It takes the form 

(j)k = Pol2k{T2,n,- ■ ■r2k+2) (21) 

and depends on 2k of the r variables only, T2,3,4,...,2fc,2fc+2- In general, only when 2k + 2 > n 
does the kth correction begin to depend on all n of the r variables. Hence the first corrections 
(which are important in practice) contain very few r's independently of n. 

The first three energy corrections €1^2,3 have a quite non-trivial property - they vanish 
at non-physical values of n = 0, 1, — It seems quite natural to conjecture that the 
correction of arbitrary order will continue to have this property so that 

where ek{n, v) is a polynomial in ra, v. Most likely there exist some physical reasons behind 
of this property, but the present author is not aware of them. 

It is worth mentioning that there is no doubt that the present perturbation theory (fT^ 
is divergent. The coefficients should grow factorially with k. However, it is not clear how 
to calculate the index of divergence. 

If (7 = in ((2|), the singular term in the potential disappears and the formulae for 
corrections simplify. This happens when i/ = 0, 1 (see (^) for which 



at 1/ = 



?)n{n — 1 



62 



2(n + 6)cj2 ' 
?,n{n - l)(25n + 6) 



at u = 1 



£2 



8(n + 6)2cj5 

9n(n - 1) (215^2 + lUn + 36) 
16(n + 6)3cu8 

^71^ {n^ — 1) 
" 4(n + 6)cj2 ' 

n2(n2 - l)(101n2 + 36) 



£3 



16(n + QYuj^ ' 

n^{n^ - l)(1026n^ + lOSSn^ + 324) 
16(n+6)3a;8 
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It is interesting to point out that at z/ = 1 the numerators of 61,2,3 depend on n^. This could 
be a general feature hold for arbitrary correction. 

For the two-body case, n = 2, the problem is reduced to a standard one-dimensional 
anharmonic oscillator where z/ = and v = 1 cases correspond to the ground state and the 
first excited state, respectively. Explicitly the corrections are. 



at = 



_ 3 _ _ 21 _ 333 



at z/ = 1, 

15 165 3915 



^1 — TT^ 1 ^2 — — „„ ^ , £3 



in agreement with the results of the calculation carried out in and Q. 

2. First excited state 

The first excited state of the gauge- rotated Calogero Hamiltonian /icai (see is char- 
acterized by, 

{n-l){l + vn) , , 

A simple analysis of equation (fTHI) shows that the first eigenfunction correction 0i G P3, 
hence it should be a third degree polynomial in r's. After substitution of such an Anzatz 
into (fTB|) simple calculations give 

n[6 + z/(5n — 6)][n + 11 + vn{n — 1)] , . 

= 4(n + 6K ' ^24) 

n n[n — 9 + z/n(n — 1)1 1 o 

~ T4.T2 H ^7 , o ^4 - — 



2(n + 6)cj 8(n + 6)cj2 4cu 

ra^ - 27ra - 54 + unin^ - 15n + 1 



16(n + 6)tu2 



■ To 



n[6 -|- z/(5n — 6)][n -|- 11 + uni^n — 1)] 
16(n + 6)u!^ 

It is worth mentioning that the correction 0i depends on two r variables only, T2^a (cf. ()18j) ). 

A similar analysis of equation (fTBjl shows that the second eigenfunction correction 02 G 
P5, hence it should be a fifth degree polynomial in r's. After substitution of such an Anzatz 



11 



into ()16|) simple calculations give the following results for the second correction 
€2 = — — ^ , r - 3n^ + 342^3 + 6471^2 + 5574n - 360 

- 2z/n(3n^ - 437n^ - 456471^ + 3666n + 2916) 



— u n 



^-^{3n^ - 73An^ - 2421n^ + 9092n - 6324) 



+ 2iy^n%n - l)(101n2 - 245n + 150) 

We will not present the explicit form of 02 due to its complexity. It is worth mentioning that 
the correction 02 depends on four r variables only, T2,3,4,6 as happens for the ground state 
(see fjl9p ). Neither ei nor e2 vanish simultaneously for some value of n. It can be shown 
that the A;th eigenfunction correction has the property (pk £ 'P2k+i- 

If = in (j2I), the singular term in the potential disappears and the formulae for 
corrections simplify. This happens when u = 0,1 (see Q) for which 

• at z/ = , (po = T2 + {n — 1)/4:U , 

_ 3n{n+ 11) 
~ 2(n + 6)cu2 ' 

3{n^ - lUn^ - 2157^2 - 1858n + 120) 

^2 — 



32(n + 6)2a;5 
at u = 1 , 00 = T"2 + (^^ - l)/4u; , 

5n2(n2 + 11) 

ei 



4(n + 6)u;2 ' 

199n6 + 36n^ + 4082^^ + 78^3 + 5463^2 - 258n - 360 



£2 



32(n + 6)2cj5 

For the two-body case, n = 2 the problem is reduced to a standard one-dimensional 
anharmonic oscillator where u = and u = 1 cases correspond to the second and third 
excited states, respectively. Exphcitly the corrections are, 

• at 1/ = 0, 

39 615 



^1 ~ ' ^2 



32uj^ 



at z/ = 1, 

75 1575 



^1 — TT^ 5 ^2 



n n 

in agreement with the results of the calculation carried out in ^] and |10[ • 
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C. Correlation functions and perturbation theory 



By purely algebraic means we calculated the first correction to the ground state energy 
(jl8|) . Making a comparison of this result with a formula for the first energy correction in 
the Rayleigh-Schroedinger perturbation theory we find ^ that, in fact, we have calculated 
the expectation value 

(0|i;p|0) _ n{n - 1)(1 + un)[6 - z/(6 - 5n)] 
~ (0|0) ~ 4(n + 6)cu2 ■ 

This expectation value is a rational function of the parameters uj^n^u 

/n n 



I \{\y^-y^'''e-^^^^y■^-y^^^'c^-'y 

n{n — 1)(1 + z/n)[6 — z/(6 — 5n)] 



4(n + 6)cu2 

where — ?/„ = Xlj"=a^ Vi ^^e domain of integration is the principal Weyl chamber. It is 
quite amazing that although each integral is a complicated combination of Euler F-functions, 
their ratio reduces to the rational function. This is a general property which appears in 
algebraic perturbation theory ^j. 



D. Variational study 

We consider a strong coupling limit A ^ oo in ((2)), which is equivalent to putting uj = 
g = 0. Following the recipe for choice of the trial functions (see e.g. [Kl), the simplest trial 
function for the ground state can be written as 

^^^^^^^ = e-(-r2)-§/3(-r2)3/2-^(a2+r|)l/2_5(.2(^.4_,|))l/2 ^ ^25) 

where a,f3,'y,S and a are variational parameters. From dimensional arguments it seems 
clear that the ground state energy should be of the form 

E„ = /(n)A^ . (26) 

For two- and three-body cases the result of calculations is 

/(2) = 0.53042 (a = 0.837, P = 0.837, a = 7 = 5 = 0) , (27) 
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/(3) = 1.17273 (a = 0.914, (3 = 0.845, a = 7 = 5 = 0) . (28) 

For the two-body case one can make a comparison with the best numerical studies, 
f{'^)numerics = 0.530362... (see e.g. Hence the simple trial function reproduces four 

significant digits in the energy. 



III. ANHARMONIC WOLFES MODEL 

Let us introduce the Hamiltonian which describes a system of three identical particles 
with two- and three-body interactions 



He 



1 ^ 



k=l 



dxl 



I 2 2 



1 



k<l 



(29) 



where u is the parameter, g = z/(z/ — 1) > — ^ and gi = fi{fi — 1) > — ^ are the coupling 
constants associated with the two-body and three-body interactions, respectively, A > is 
an anharmonic coupling constant and the factor 1/36 is introduced for convenience. We call 
this system the G2 anharmonic oscillator. 

At A = the Hamiltonian ()29j) becomes the Hamiltonian of the rational G2 model 



which was introduced for the first time by Wolfes 



Reduction method 
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and later obtained in the Hamiltonian 



. Its ground state is given by 



¥^\x) = (Af)(x))^(Af^(x))'^e-5-SLi-I , Eo = -uj{l + 2u + 2/i) , (30) 



where A^{\x) = ni<j \xi - Xj\ and A^^'i^) = Ilt<j, ij^k + ~ "^^kl- 

An interesting observation is that all fourth order permutationally symmetric and 
translation invariant polynomials correspond to two body interactions because 

{Xi + X2- 2X3)^ + {Xi + Xs- 2X2Y + {X2 + X3 - 2xiY = 

= 9 [(xi - X2)' + (xi - xs)^ + ix2 - xa)^] , 



(r), 



and 



(xi - X2)^(xi - xa)^ + (xi - X2)^(x2 - xa)^ + (xi - Xa)^(x2 - xa)^ 
= 1/2 [(xi - X2)^ + (xi - xa)^ + (X2 - xa)^ . 
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This leads to the important conclusion that a general fourth degree polynomial transla- 
tion invariant potential reduces to two body interactions. Therefore, the Hamiltonian ()29|) 
describes the most general permutationally symmetric and translationally invariant anhar- 
monic oscillator associated with the G2 rational model with fourth order polynomial anhar- 
monicity. 

Let us make a gauge rotation of the Hamiltonian ()29p with the ground state eigenfunction 
(EOl), 

he, = 2i¥;\x))-\nG. - Eo)¥;\x) . (31) 

The result can be written in terms of two relative coordinates and the center-of-mass coor- 
dinate X. 

Now let us take the Perelomov relative coordinates (0) and introduce new permutationally 
symmetric relative coordinates, 

Ai = yl + yl + yl , A2 = yMyl , (32) 
with the condition —y^ = Z/i + 1/2 (cf- 0)- Making in a change of variables 

(a;i,X2,X3) (F, Ai, A2) 

and separating the center-of-mass motion (and then omitting it), the remaining part of the 
Hamiltonian pip takes the form 



hG. = -^Xidl,, - 24X^81,,^ - I8XIX2& 



A2 A2 



+ {4uXi - 4[1 + 3(^ + u)]}dx, + [l2uX2 - 9(1 + 2z/)A?] dx, + XXj . (33) 

This is the algebraic form of the G2 anharmonic model (cf. © at = 3). This Hamiltonian 
possesses a remarkable property - among eigenfunctions there exists a family which depends 
on the variable Ai only (!). The ground state belongs to this family. In order to find the 
eigenfunctions depending on Ai only it is necessary to solve the spectral problem for the 
operator 

he = -4X^^l^^ + {AuXi - 4[1 + 3(/i + u)]} d^, + AA? . (34) 
By making a gauge rotation the operator (jH^ can be reduced to the two-body Hamiltonian 



^--2 2^d^ + Y^'''''''> + (x, - X2y +8 ^"^^ ' 

i=l ^ 



\4 



y4i— rational model 
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Similarly to what was done for the An anharmonic many-body oscillator in Section 1.2, 
one can develop perturbation theory in powers of A for the Hamiltonian (|^^ taken in the 
algebraic form 

IV. CONCLUSION 

We introduced an anharmonic perturbation of two completely-integrable and exactly- 
solvable systems, which are in fact anharmonic many-body oscillators. It is not clear 
that these systems remain integrable or whether the anharmonic terms break this feature. 
However, the calculation of perturbation corrections is not influenced by existence or non- 
existence of integrability. Perhaps, it is interesting for the A^-anharmonic oscillator to study 
the limit n oo and a field-theoretic limit. Another interesting question is about the exis- 
tence of the quasi-exact ly- solvable anharmonic generalizations of the Calogero and Wolves 
models other than those found in jl^l . 
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